Hitting time for the continuous quantum walk 
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We define the hitting (or absorbing) time for the case of continuous quantum walks by measuring 
the walk at random times, according to a Poisson process with measurement rate A. From this 
definition we derive an explicit formula for the hitting time, and explore its dependence on the 
measurement rate. As the measurement rate goes to either or infinity the hitting time diverges; 
the first divergence reflects the weakness of the measurement, while the second limit results from 
the Quantum Zeno effect. Continuous-time quantum walks, like discrete-time quantum walks but 
unlike classical random walks, can have infinite hitting times. We present several conditions for 
existence of infinite hitting times, and discuss the connection between infinite hitting times and 
graph symmetry. 



I. INTRODUCTION 



There are two main types of quantum walks: 
continuous-time and discrete-time quantum walks. 
Discrete-time quantum walks evolve by the application 
of a unitary evolution operator at discrete time inter- 
vals, and continuous-time walks evolve under a (usually 
time-independent) Hamiltonian. Continuous-time quan- 
tum walks have been defined by Farhi and Gutmann in [l[ 
as a quantized version of continuous-time classical ran- 
dom walks. Classical random walks are used in computer 
science to design probabilistic algorithms for computa- 
tional problems most notably for 3-satisfiability (3-SAT) 
Q . In a similar vein, quantum walks provide a framework 
for the design of quantum algorithms. As such quantum 
walks have been used in many quantum algorithms such 
as element distinctness Jj|, matrix product verification 
triangle finding Q and group commutativity test- 
ing [||. Recently, a quantum algorithm for evaluating 
NAND trees has been proposed which uses a quantum 
walk as a part of the algorithm In order to be able 
to understand how to better use quantum walks for algo- 
rithms, we need to study the properties of these walks. 
There have been many papers which study the behav- 
ior of quantum walks for various graphs. For example 
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quantum walks on the line have been examined for the 
continuous-time case in Refs. HHLIsI an d f° r the discrete- 
time case in [hJ [H d GJ, [OfT The iV-cycle i s treated in 
US [3, and the hypercube in [H 03, Gj, SS EU . Quan- 
tum walks on general undirected grap hs are defined in 
[HI, HH, and on directed graphs in |24j . Kendon [25[ has 
a recent review of the work done in this field so far, fo- 
cusing mainly on decoherence. Other reviews include an 
introductory review by Kempe in [2(| , and a review from 
the perspective of algorithms by Ambainis in [23| . 

Different quantitative characterizations of quantum 
walks are defined by analogy to the classical ones, such 
as mixing times, hitting (absorbing times), correlation 
times, etc. Often for this purpose the evolution of 

the quantum walk must be modified, to include not only 
the Hamiltonian evolution, but also a measurement pro- 
cess to extract information about the current state of the 
walk. There is a natural way to introduce such a mea- 
surement process in the discrete case: namely, a mea- 
surement is made after each step of unitary evolution. 
The outcome from this measurement is used in defining 
the characteristic time scale in question. In the case of 
the continuous-time walk, such a natural definition of a 
measured walk does not exist. There is no intrinsic time 
step after which we can perform the measurement. Clas- 
sically this is not a difficulty, because measurements do 
not disturb the state of the system. The quantum case 
is quite different. If we choose the measurement times 
arbitrarily, they can either be too long or too short with 
respect to the unitary evolution of the quantum walk. 
We can either miss important details in the evolution by 
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measuring too infrequently, or overly distort the unitary 
evolution by measuring too often. In the limiting case, 
we can completely freeze the evolution by the Quantum 
Zcno effect |27j |. 

Hitting times for discrete-time quan tum walks have 
been defined and analyzed in [l9l. 1 20t| - The effect on 
mixing times of making random measurements, and its 
possible algorithmic applications, has been studied in 
[28l. [2^. f30(| . In this paper, we introduce a measurement 
process for the continuous-time quantum walk which 
gives rise to a definition and an analytical formula for 
the hitting time as function of the measurement rate 
(or equivalently, measurement strength) . We explore the 
limits of measuring too weakly or too strongly, and show 
that the hitting time diverges in either case. This sug- 
gests the existence of an optimum rate of measurement, 
which depends on the unitary dynamics of the particular 
walk. 

We also show another difference from hitting times for 
classical random walks. In the classical case, a random 
walk on a finite connected graph always leads to finite 
hitting time for any vertex. This is not true for the 
quantum case. The existence of infinite hitting times 
has been argued for discrete-time and continuous-time 
quantum walks in [20| ; in this paper we show this explic- 
itly for the continuous-time quantum walks based on the 
definition of hitting time that we give, and derive condi- 
tions for the existence of infinite hitting times. Another 
sufficient condition that we prove is that if the comple- 
mentary graph is not connected, this automatically leads 
to infinite hitting times for the continuous-time quantum 
walk on the original graph. 

The paper is organized as follows. In section II, we 
describe how to introduce the measurement process, and 
derive formulas for the hitting time and probability. In 
section III, we give a condition for existence of infinite 
hitting times and prove that the hitting time diverges 
when the measurement rate goes to zero or infinity. In 
section IV, we give examples for the hitting times for 
certain graphs as a function of the measurement rate. In 
section V, we give another sufficient condition for infinite 
hitting times. A discussion follows in section VI. 



II. HITTING TIME DEFINITION FOR THE 
CONTINUOUS QUANTUM WALK 

We want to define hitting time for the continuous uni- 
tary evolution on undirected graph r(V, E), where V is 
the set of vertices and E is the set of edges. Two ver- 
tices v\,V2 € V are connected if there exists an edge 
e = {^1,^2} € E (here {^1,^2} should be taken as the 
unordered pair or set of the two vertices V\ and ^2). 
Corresponding to the graph T, we assign the Hilbcrt 
space TLt = l 2 (V). The vertex states in that Hilbcrt 
space are just labeled by the vertices of the graph — for 
v € V, \v) £ H.. They form an orthonormal basis for Tlr'- 

(^nl^m) — 3 nm . 



The hitting time for classical random walks is defined 
naturally as the average time to find the walk in a spe- 
cific vertex. When we turn to the quantum case the walk 
on the graph is not defined as a stochastic process on the 
vertices of the graph but as the unitary evolution of a 
closed quantum system with a Hilbcrt space defined as 
above. In order to be able to tell when the quantum walk 
has reached a vertex, we need to measure the system in 
order to gain information about the current state of the 
system. There are several reasonable ways to do that 
in the continuous-time case. We could perform strong 
measurements periodically with some fixed but arbitrary 
period T. This is not unlike in the discrete case, in which 
the period T is given naturally by the walk itself: a mea- 
surement is performed after each unitary evolution step. 
This way to perform a measurement in our case is unsat- 
isfactory. We have no way to know how to choose T. If 
we choose it too small we could introduce too much deco- 
herence, effectively masking the unitary evolution of the 
walk, or even worse, freezing it. If T is too large then we 
can miss the moment when the walk actually reaches the 
final vertex. And in general, the unitary transformation 
between measurements can be complicated and difficult 
to work out (unlike the discrete-time case) . 

Another way to measure the system is through strong 
measurements but performed at random times. The mea- 
surement times arc chosen according to some probability 
distribution with some measurement rate. The advan- 
tage is that we don't introduce an artificial periodicity 
into the dynamics, and it allows one to calculate aver- 
aged effects over different measurement patterns. The 
disadvantage is that it is still necessary to introduce a 
time scale for the measurements, this time given by the 
rate at which measurements are performed. 

A third way to measure the system is using "weak" 
measurements, analogous to the way photodetection is 
described. In a small time period 5t we perform a mea- 
surement which either allows the system to evolve unitar- 
ily with a probability 1 — e, or performs a measurement 
to determine whether it has reached the final state with 
a probability e. In this case, the evolution is unitary for 
most of the time, with jumps at random time when the 
measurement is performed. This case and the previous 
one are actually equivalent — the values of St and e deter- 
mine the measurement rate A — but they give a somewhat 
different intuition about how to look at the measurement 
procedure. In the second case, it is natural to describe 
the random times is by a Poisson process with a given 
rate. In the third case the role of a rate is played by 
the strength of the "weak" measurement. In his work on 
mixing times [3(|, Peter Richter argued that the qual- 
itative and even quantitative behavior of the system is 
not too sensitive to the exact details of the measurement 
scheme. This suggests that the first choice above might 
be as good as the the other two, but it still introduces a 
discrete structure which is not desirable in dealing with 
continuous evolution. In the following, we explore a mea- 
surement scheme described in the second and third cases. 
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We do a measurement to check whether the system is 
in the final state \vf) , given by the measurement oper- 
ators {Pf,Q f } where P f = \v f ){vf\, Q f = I - P f . We 
have to specify the times when we perform the measure- 
ments. We will measure the system at random times 
distributed according to a Poisson process X t with rate 
A > 0. Each time we observe a jump in the Poisson pro- 
cess we measure the system. Between the moments at 
which wc perform the measurements the system evolves 
unitarily with a Hamiltonian H = —"/L, where L is the 
discrete version of the continuous Laplacian V 2 . In our 
case it is given by L = A — D, where D is a diagonal 
matrix in the basis spanned by the vertex states with the 
degree of each vertex along the diagonal, and A is the 
adjacency matrix of the graph [l], [31(. In this paper we 
take 7 = 1. If the degree of the vertex v n is d n then we 
have the following representation of D and A 



where 



D = y^d n \v n )(v n \ 

n 

A = ^a„ m \v n )(v m \ 



1, if {v„,v m } G E, 
0, if {v n ,v m } g E. 



(1) 

(2) 

(3) 



is the adjacency matrix of the graph T(V,E). 

Let to = (ti, ■■•) be a sequence of random times when 
the jumps of the Poisson process are observed, with (t n G 
R, < ti < t 2 < ■■■)■ For convenience we will take to = 0. 
The sequences u> belong to a probability space (fl, P) 
on which the Poisson process Xt is defined: 



I:Rxfl^ {0,1,2,...}, 

X t (uj) =ra, if t G [t n ,t n+1 ). 



(4) 



Here is the set of all sequences of random times to, $ 
is the cr-algebra, generated by the Poisson process. The 
probability measure P on £1 is the one induced by the 
Poisson process. (For reference, see (3"3.l33T|.) 

For each sequence u> G we define the hitting time as 



We would like to find an analytical formula for t),. 
From the definition of p n , we have 



n-l 



m— 1 
n-l 



On = Tr I P f J] U-Kt m+1 -t m ) HQ \ e -it lH p .^H x 



n(>< 



i(t m+1 -t m )H' 



(7) 



m— 1 



The arrow above the products signify whether the oper- 
ators entering the products are ordered from left to right 

< 

or vice versa, in other words [Tjt/i = U n U n _\...Ui and 

fcZut = u l u 2 ...u n . 

We now introduce superoperators Ui and Q f . defined 
by 



-m Xe itH 



U- t (X) = e 
Qf(X) = Q f XQf 



(8) 
(9) 



and use them to rewrite ([7|) as 

Pn = Tr {PfU~ tn oQfO Z4 n l o Q f o ... o U h ( Pl )} . (10) 

We want to express the sum, |5}, as a function of the 
{t„}. We do that by adding, subtracting and rearranging 
terms in the sum, assuming that it is absolutely conver- 
gent: 

Tu, = hpi + t 2 P2 + hPz + UP4 + ... 

= hpi + t 2 P2 - hP2 + hP2 + t 3 p 3 - t 2 P3 

+ tiPz - hpz + hp 3 + tiPi - ... 

= hiPl +P2+P3 + ••■) + (*2 - h)(P2 +P3 + ■••) 
+ (t 3 -t 2 )(P3+Pi + ...) + .• ■ 



k—1 ri—k 



(11) 



Because the {t n } are the event times of a Poisson process, 
the interval times {i n } are independent and identically 
random variables, exponentially distributed with param- 
eter A and a probability density function given by 



n=l 



tnPn 



(5) 



where p n is the probability to find the system in the fi- 
nal state at time t n given that the system wasn't mea- 
sured to be in the final state in any of the previous times 
t n -\,...,t\. (For reference, see [IS HH-) Define the in- 
tervals between jumps tj-i and tj as ij = tj —tj-i. We 
want to average r w over all possible trajectories u of the 
Poisson process, and take this as our definition for the 
hitting time: 



T h = E P (r w ) = / r w dP(w) 



(6) 



JtnK ' [0, t<0. 
Knowing that, we reexpress formula ([6]): 



(oo 00 
n / ^ 
tin 



-At, 



Then 



Th 



OO OO / oo „ \ 

ee n w 

7. 1 m 7„ \ 1 1 J I 



(12) 



(13) 



(14) 



k—1 n—k \ l — l 



4 



In the above expression there are two types of integrals: 



A{X) = j dt\e- M Ut{X), 
o 

oo 

B(X) = J dt\e- xi UA- t {X). 
o 



(15) 



(16) 



Integrating by parts we get the following equations for 
the operators A and B 



A + j[H,A] = X, 

B + t [HjB ] = ~A(X), 



(17) 
(18) 



where A(X) is the solution to the first equation. (We 
will prove that this solution exists below.) Defining the 
superoperator 



C X (X)=X + -[H,X], 



we rewrite these equations as 



C X {A)=X, 
£ X (B) = \- 1 A(X). 



(19) 
(20) 



We want to prove that the superoperator C\ is invertiblc 
when A is a real number. For this we need to know how 
the adjoint of a superoperator is defined with respect to 
the Hilbcrt-Schmidt inner product for operators, 

(X,Y} HS = Ti(X^Y). (21) 

Using this inner product, we see that if 

C(X) = Y. c " C " XD l 

n 

the adjoint of C(X) is given by 

cHx) = Y / <cixD n . 

n 

From that it follows that C is a normal superoperator: 



4 o C x - C x o 4 = 0. 



(22) 



This means that C\ is diagonalizable. If X n is an eigen- 
vector of C\, then X„ is an eigenvector of the hermi- 
tian and anti-hermitian parts of C\ separately, which are 
given by 

C^{X)= l -{Cx+C\){X)=l{X)=X 



and 



Ci{X) = -{C x -C\){X) = \[H,X], 



respectively. Let us denote the eigenvalue of 4 cor ~ 
responding to X n by ix n [x n £ R because >4 is anti- 
hermitian). Then C x (X n ) = (1 + ix n )X n ^ 0. This 
proves that each eigenvalue of C\ is nonzero, and that 
C\ is invertible. 

The solutions to equations (fTT)) and (TT5]) are 



A = £-\X), 
B = \- 1 C- X \X). 

Substituting these in pT]) we get 



(23) 
(24) 



OO OO 



k—1 n—k 

°{Qf°£-x l ) k ~ l {Pi)} 

OO OO 

1 n—k 

o£- 1 o(Q / o£- 1 ) fe - 1 (p I )} 
£f;A- 1 Tr{p / £^o(Q / o£- 1 ) i c 



fc=i \ 



k=0 1=0 



o£- 1 o(Q / o£- 1 ) fe ( Pl )} 



(25) 



Note that the eigenvalues of C\ are always greater than 
or equal to 1 in absolute value, from which it follows that 
the eigenvalues of Qf o C^ 1 are all less than or equal to 
1 in absolute value. If all the eigenvalues are strictly less 
than 1 in absolute value, the following sum exists: 

OO 

Substituting this into (f25j) and denoting J\f\ = C\ — Qf, 
we get the following formula for the hitting time: 



r h = A _1 Tr {PfM^ 2 {pi)} 



(26) 



This formula is closely analogous to the formula for the 
hitting time derived in [2(| Hi} for the case of a discrete- 
time quantum walk. If Qf o C^ 1 has any eigenvalues 
equal to 1, the inverse in the above formula should be 
thought of as a pseudoinverse. 

Another quantity that may be defined is the total prob- 
ability to ever hit the final vertex [2l|, [34| : 



Ph = ^2pn- (27) 

n=l 

We can derive a formula similar to (|26[) for ph : 

Ph = TT{Pf(C x -Qf)- 1 (p l )}. (28) 

When the hitting time is not infinite, or equivalcntly 
when the superoperator C\ — Qf is invertiblc, ph = 1- 
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In the case of infinite hitting time, we replace the inverse 
with a pseudoinverse, as before. 

There is another way to derive formula (f2B|) : by look- 
ing at this procedure as an iterated weak measurement 
(case 3 that we discussed above). Instead of summing 
over all trajectories of the Poisson process at each time 
period St, we perform a generalized measurement with 
measurement operators 

M = ^/l-e 2 e- lStH , 

Ah = sP f , (29) 

M 2 = sQ f . 

These operators form a complete measurement as 
Y^i=o M^Mi = I. The measurement is weak (in the par- 
ticular sense of giving little information about the sys- 
tem on average) when e <C 1. Let us define a positive 
matrix p c describing the state of the system at time t 
conditioned on the assumption that outcome "2" has not 
occurred up to this time. We measure the system repeat- 
edly at intervals of time St, using the same measurement 
operators, and if we don't observe outcome "2" the state 
of the system is described by the following matrix: 

l 

P c {t + 5t) = Y J M lP c (t)hll (30) 

i=0 

We expand in powers of the small parameter e and take 
the limit e — > and St — ► 0, keeping the ratio e/VSt con- 
stant, and obtain a master equation for p c (t). This gives 
the connection between the strength of the measurement 
e and the measurement rate A = lim^^o e 2 /St. After we 
expand to second order of e we get 

p c (t + St) =(1 - e 2 )(l - iStH)p c {\ + iStH) 
+ e 2 Q f p c Qf + 0(e 3 ) 
=p c (t)-iSt[H, p c ] 

-e 2 (p c (t)-Q fP c Q f ) + 0(e 3 ). (31) 

Taking the limit St ^ and using the fact that 

Um p c (t + St)-p%t) = dp? 
<5t->o St dt ' 

e 2 

lim — = A 

<5i-»0 St 

we arrive at the master equation for p c : 
do c 

= -i[H, P c ] -X(p c - Q f p c Q f ) = -W X ( P C ). (32) 

Note that p c is positive, but not normalized; the trace 
of p c is the probability that measurement result "2" has 
not been seen up until time t. The total probability to 
hit the final vertex ph and the hitting time Th are given 



by: 

OO 

Ph = \ J Tr{P f p c (t)}dt, (33) 
o 

OO 

r h = \J m{p fP c (t)}dt. (34) 



Substituting the solution of (|32p 

p c (0= e - At ^(p c (0)). 

in ([33| and (|34l) and integrating by parts we obtain for- 
mulas (|26p and (|28p for the total probability to hit and 
the hitting time. 

III. CONDITIONS FOR EXISTENCE OF 
INFINITE HITTING TIMES 

We want to prove that the existence of infinite hitting 
times is equivalent to the non-invertibility of the super- 
operator C\ — Qf. We will need the following definitions 

[Mill. 

Definition 1. A matrix pencil A + sB (where A and B 
are n x n matrices and s is a complex number) is said to 
be regular if there exists at least one complex s for which 
the pencil is nonsingular. 

Definition 2. A complex number s is a finite eigenvalue 
of the regular matrix pencil A + sB if det(A + sB) = 0. 

Definition 3. The regular matrix pencil A + sB is said 
to have an infinite eigenvalue if B is a singular matrix. 

In discrete-time quantum walks, infinite hitting times 
were observed for specific graphs and initial states p51 . 
[2l| . This occurs when, starting from the initial state, 
the total probability to ever find the walk at the final 
vertex is less than 1. It has been argued in the pa- 
pers above that infinite hitting times occur given that 
the graph has a symmetry group that leads to a degen- 
erate Hamiltonian of the quantum walk. The symmetry 
of the group leads to splitting the Hilbert space into in- 
variant subspaces under the action of the group, on each 
of which the group acts with one of its irreducible repre- 
sentations. The evolutionary operator for the quantum 
walk leaves these subspace invariant under its action be- 
cause the symmetry group of the graph is necessarily the 
symmetry group for the Hamiltonian. Thus if the final 
state is in one invariant subspace, but the initial state 
does not lie entirely in the same subspace, there will be 
a nonzero probability to never hit the final state. For 
such a situation to occur the final state, which we always 
assume to be a vertex state, must lie entirely in one of 
those invariant subspace. A sufficient condition for that 
in the case of a discrete-time quantum walk on a regular 
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graph is the presence of an irreducible representation of 
the symmetry group in the Hamiltonian with a dimen- 
sion larger than the dimension of the coin space. For the 
continuous-time quantum walk the equivalent condition 
is the presence of an irreducible representation with di- 
mension larger than one. As Abelian groups always have 
one-dimensional representations, we might expect that 
in order to have infinite hitting times we need to have a 
symmetry group that is not Abelian. We will show that 
this is not true. The symmetry of the graph can lead 
to infinite hitting times even if the Hamiltonian is not 
degenerate, as is the case when the symmetry group of 
the graph is Abelian. Having a non- Abelian symmetry 
group is a sufficient, but not a necessary, condition. 

Another condition for the existence of infinite hitting 
times is connected to the invcrtibility of the supcropcra- 
tor A/a which enters formulas ([26]) and ([28]) . Consider all 
operators X such that [H, X] = and PfX = XP f = 0, 
and denote the projector on the linear subspace of all 
such operators by V . We will prove that V ^ if and 
only if N\ = C\ — Qf is not regular. 

HV^O then choose X such that V(X) = X ^ 0. 
Then [H, X] = and from P f X = XP f = follows that 
Qf(X) = X. Thus Af\(X) = which means that N\ is 
singular for every A and thus not regular. 

If M\ is not regular then it is non-invertible for any A. 
Let us fix A to be real and different from 0. There exist 
al/O such that N\(X) = 0. We have already proven 
that C\ is invertiblc for real A. Then 

Mx(X)=Cxo(l-C^oQ f )(X)^0 (35) 

and thus 

C~ x 1 oQ f {X)=X. (36) 

Taking into account that all eigenvalues of C\ are greater 
or equal to 1 in absolute value the above equality is true 
only if 

CT X \X)=X, (37) 
Q f (X) = X, (38) 
which are equivalent to 

[H,X]=0, (39) 
P f X = XP f = 0. (40) 

This means that V is nonzero. 

Let's explicitly give the form of the projective super- 
operator V in terms of the projectors on the eigenspaces 
of the Hamiltonian H. For this purpose we have to de- 
fine the intersection operation n for orthogonal projec- 
tions. For any two orthogonal projection operators P\ 
and P2, Pi (~1 P2 will denote the orthogonal projector on 
the subspace which is the intersection of the subspaces 
onto which Pi and P2 project. Let the Hamiltonian H 
have the following decomposition: 

r 
i=l 



where Ei are the eigenvalues of H (£", ^ Ej for i 7^ j), 
Pi are the projectors onto the eigenspaces of H cor- 
responding to eigenvalues Ei. Since H is Hermitian, 
PiPj = SijPi, and Tr Pi = di is the multiplicity of the 
Ei eigenvalue. The projector V is then 

r 

v(x) = J2(P*riQf)x(p l nQ f ). (4i) 

From this form it is easy to see that V is a completely 
positive superoperator. As such, if it is different from 
0, then there must exist a density matrix p such that 
■p(p) = p. If the walk begins in such a state p, it will 
never arrive at the final vertex. 

Now we will prove that if Af\ is regular as a matrix 
pencil then all its eigenvalues lie on the imaginary axis. 
Let's assume that Af\ is invertible for some A ^ and 
Af\ is non-invertible for some Ao 7^ 0, Ao 7^ A. Then 
N\{X) ^ for all X ^ and there exists X ^ such 
that Af\ (X ) = 0. Then (Af\ - AA Ao )(A ) = i(l/X - 
1/X )[H,X ] ^ and therefore [H, X Q ] ^ 0. Analogously 
(W A -AoAA Ao )(A ) = (A-A )(I-Q/)(X ) + therefore 
{I — Qf)(X ) 7^ 0. As 1 — Qf is a projector, it follows 
that 

(X ,(T-Q f ){X )} HS ^0. (42) 

Taking into account that 1 — Qf and H(-) = [H,-] are 
both Hermitian supcropcrators, (X , (1 — Qf)(X ))HS 
and (X , H(X ))hs a re both real numbers. Denoting 
Aq = i?e(l/Ao) and X l = im(l/Ao) we have 

(X ,(J- Qf-\iH)(X )) HS + i\ r (X ,H(X )) HS = 0. 

. (43) 

This equality is only possible if both the real and imag- 
inary parts vanish, implying that Aq = and hence 
i?e(A ) = 0. 

Using this result, we will now prove that if J\f\ is a regu- 
lar matrix pencil, and thus doesn't have infinite eigenval- 
ues, then the hitting time behaves regularly as a func- 
tion of A on the real line: it won't diverge for any real A 
except when A goes to or infinity. Physically, this means 
that when we measure either very weakly or very strongly 
we never find the particle in the final vertex. The first 
limit is easy to understand — if we never measure, we will 
never find the particle anywhere. The second limit — A 
going to infinity — corresponds to the Quantum Zeno ef- 
fect, in which the evolution of the system is restricted to 
a subspace orthogonal to the final vertex. 

To prove this conclusion, we represent supcropcrators 
as matrices using the following isomorphism: 

0:C(-)= $3 c„G n (■)£>£ 

n 

— »^(C) = C = J2 C nCn ® D* n . (44) 



7 



Now we represent the superoperator pencil j\f\ by the 
matrix pencil 



= I®I-Q f ®Q)-\{H®I-I®H*). (45) 



By assumption this matrix pencil is regular. Every reg- 
ular matrix pencil A + sB has the following canonical 
form: 

A + pB = T{N {mi \...,N {m *\ 

J^(p ni ),...,J^(p nq )}S, (46) 

where T and S are invertible matrices, constant with 
respect to p, and {7V( mi \ is a 

block-diagonal matrix where the blocks and J^ n > 

are square matrices of order m and n respectively of the 
form 



(1 ix 
1 p 
1 




\o 



\ 





1 p 

o I J 



1 



( M + Mi 

/x + /xj 






1 

M + Mi 






(47) 



\ 






n + m i 

p + pi ) 

(48) 



or more succinctly = 7 (fc ) + pK^ and = 

(p + po )I^- n ' +K^ n \ where /' m - ) is the the identity matrix 
of order m and K^ m ^ is a m x m matrix with "Is" imme- 
diately above the diagonal and "0s" everywhere else. The 
]\[( m ) blocks are present when the matrix pencil has infi- 
nite eigenvalues and the j( n ' blocks correspond to finite 
eigenvalues. 

We want to examine the behavior of the inverse of a 
regular matrix pencil when p approaches one of its eigen- 
values. Assume that p$ is a finite eigenvalue of the pencil 
and the corresponding n x n block to that eigenvalue is 
J(po). The inverse of this block is given by 



3=0 



j+1 



(49) 



In the above K° = I and K = ffW . The inverses of all 
blocks that do not correspond to the eigenvalue po will 
have regular behavior when p approaches — hq. We want 
to examine the behavior of when p goes to infinity. 
Analogously the inverse of this block is given by 



(50) 



The inverse of the blocks corresponding to finite eigenval- 
ues will have regular behavior when p approaches infinity. 
As the matrix pencil 



N, 



Ni/ M = I®I-Q f ®Q* f + ip(H(g>I-I(g)H*) (51) 



is regular and has both finite {p = 0) and infinite eigen- 
values (because both matrices / <E) I — Q / <£> Q*f and 
i {H ® I — I (g> H*) are singular), both types of blocks 
an d J(0 are present in its normal form. If we ex- 
press formula (|26[) in terms of matrices and vectors with 
p = 1/ A we get an analogous expression 



(52) 



where PJ and p^ are the vectorized versions of the ma- 
trices Pf and pi. 

When p goes to we can see from formula (|4"9"]) with 
Po = that the asymptotic behavior of is given by 



T h =Tl(p) + -PJ -(S-'PoT 

p 



(53) 



where rj^(p) is a function which is regular in a neighbor- 
hood of p = 0, T and S are the invertible matrices in the 
canonical form (|46|) of the matrix pencil N M and Pq is the 
projector on the eigenspace with eigenvalue p = 0. Now 
as long as P v f ■ (S*" 1 PqT- 1 ) 2 p\ ^ 0, r h will go to infinity 
when p goes to 0(A going to infinity) no matter whether 



terms of higher order, O y-^j, are present or not. 

Analogously, when p goes to infinity, 77, has the asymp- 
totic behavior 



Th 



pP] ■ (S-^T 



0(p 2 ), 



(54) 



where P^ is the projector on the eigenspace with infinite 
eigenvalue. Here again as long as PJ ■ (S~ 1 P 00 T~ 1 ) 2 Pi ^ 
0, Th will go to infinity when p goes to infinity (A going 
to 0) no matter whether terms of higher order, 0(p 2 ) 1 
are present or not. 

As we shall see in the next section the hitting time in 
the examples that we will give below has the following 
form as a function of A: 



Th 



= T (1) A- 



A 



(55) 



where the constants rn) and T(_i) depend on the partic 
ular graph. 



IV. EXAMPLES 

In this section, we will consider as examples the graphs 
in Fig. [U using the labeling of the vertices given in the 
figure when necessary. We can see examples of infinite 
hitting times {jph ^ 0) for the graphs L3, A'3, L4, KL31 
and 5*4. 

We can describe the probability to hit and the hit- 
ting time Th in another way, by specifying two operators 
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FIG. 1: Graph examples (with assigned vertex labels) 



?fi\(T,Vf) and fi\(T, Vf), and calculating their expecta- 
tions in the initial state: 



p h = Tr{ s p A (r,t; / )pJ, 
T h = Tr{f)\(r,Vf)pi}. 



(56) 
(57) 



Here pi is the density matrix describing the initial state of 
the system. These equations follow from formulas (|26p 
and l|28p. respectively. By using the definition of the 
Hilbert-Schmidt inner product, we derive the following 
formulas for ^\(T,Vf) and $j\(T, vj): 



y x (T,v f )=((C x -Q f )- 1 )\p f ), 



Si x (T,Vf) = j[(Cx 



(Pf)- 



(58) 
(59) 



In the following, we will show the operators ^a(I\ vt) 
and S)\(T,Vf) in the vertex state basis for each of the 
graphs on Figure [l] and give a brief discussion of the 
quantum walk on each graph. It is useful to describe the 
hitting probability and time in terms of these matrices, 
because they give the result for any starting state. 



The dependence of the hitting time for vertex v\ can in- 
clude two terms: the 2/A term diverges as A — * 0, which 
simply represents the increasing time it takes to find the 
particle as the measurement rate goes to zero; if the sys- 
tem starts at vertex there is also a A/2 term, which 
diverges as A — > oo because of the Quantum Zeno effect: 
as the measurement rate increases, we can "freeze" the 
system's evolution. 



2. Example 2 

The situation is different in the case of the L3 graph. 
The graph again has symmetry group C2 , and the Hamil- 
tonian has no degeneracies. Despite that, however, one 
of the three energy eigenstates has zero overlap with the 
U2 vertex: (l/y/2, 0, — l/v2)- This means that even with- 
out degeneracy, there is an infinite hitting time for the 
final vertex w 2 . This is not accidental and the symmetry 
of the graph is still responsible for the existence of this 
infinite hitting time. Under the action of C2 each energy 
eigenstate |e^) will have to be either symmetric or anti- 
symmetric. The Hilbert space thus splits into a symmet- 
ric and anti-symmetric subspaces which are orthogonal 
to each other. As the vertex state \v2) is obviously sym- 
metric under the action of the group it will be orthogonal 
to the anti-symmetric subspace. This is what leads to an 
infinite hitting time. This is a general observation for 
any graph that has C'2 as a symmetry group and vertex 
state that are left invariant under the action of the group. 
There arc no infinite hitting times for reaching vertices 
i>i and V3. Wc can sec all these properties by examining 
the matrices ty\ and 




(62) 









1 


- i 
i a 


-2X + 1 

* + x 

1 i 


A 
2 


2A 
3A 


2 2 


2A 2 




2 



_1_ 

2A 



(63) 



1. Example 1 



The existence of infinite hitting time for reaching vi can 
easily be seen from the ^ matrix for V2'- 



The graph K2 has the symmetry group C2. As this 
group is Abelian, the Hamiltonian of this graph is non- 
degenerate. The two eigenvectors have nonzero overlap 
with both vertex states, and therefore there can be no 
infinite hitting times, as we can see from the matrices ?P 
and f): 



^x{K2,vx) 



1 
1 



f)x{K 2 ,v 1 )= ( \ 2 x A 

A A 



(60) 
(61) 



%x(L3,v 2 ) 
^x(Ls,v 2 ) 






8A 

■F 1 


4A 2 


4 V 


8A 








F 1 




1 4 
9 


A 

1 


■l A 1 


1 4 
9 


8 ^ 


8A 


4A 


4 8' 


8A 



(64) 
(65) 



As ^^(^3,^2) is not the identity, there must be initial 
states that will result in a less than 1 probability to hit 
w 2 - For example, this will be true for any initial state 
which is a superposition of states |i>i) and |t>3) . 
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3. Example 3 

The graph has symmetry group D 3 , and its Hamil- 
tonian is degenerate. It has infinite hitting times to hit 
any vertex. If we calculate the <p and 5} matrices for this 
graph, we discover a new property of these hitting time 
and hitting probability matrices. The graph K 3 is not 
isomorphic to L 3 , but its ^ and f) matrices for any ver- 
tex of the K 3 are also given by (|64|) and (f65j) (or their 
appropriate cyclic permutations). This is because the 
Hamiltonians of the L 3 and K 3 graphs commute. We 
will observe the same kind of behavior below for other 
graphs with commuting Hamiltonians. 



4- Example 4 

The quantum walk on the graph L4 has the same qual- 
itative behavior as the walk on L 2 . They both have the 
same symmetry group, C 2 , as does the L 3 graph. But in 
the case of L4, as in the case of L2, there are no infinite 
hitting times. 



5. Examples 5 and 6 

We will examine the graphs KLz,\ and S4 together, 
because it turns out that their behavior is closely re- 
lated. The graph KL 3< \ again has C 2 for its symmetry 
group. The Hamiltonian is nondegenerate, but there are 
infinite hitting times for the vertices v\ and v 2 , due to the 
existence of an eigenvector which vanishes on those two 
vertices: (0, 0, 1/v2j — This is quite analogous to 

the case of the graph L 3 . In general, graphs with C 2 sym- 
metry will have infinite hitting times for hitting vertices 
that are fixed points under the action of the symmetry 
group. 

The graph S4 has D 3 as a symmetry group. Its Hamil- 
tonian is degenerate, and it has infinite hitting times to 
hit any vertex. It turns out that the matrices and 
f) for hitting vertices v\ and v 2 coincide with the same 
matrices for the graph KL 3> \; 



y x (KL 3A ,v 2 )=y x (S A ,v 2 ) 

(\ H\ 

10 

ion 

1 nil 

\ 3 U 3 3 

fi x (KL 3>1 , v 2 ) = fi\(S 4 , v 2 ) 



(68) 



/ 



_8_ J_ 

9A 3A 



_8_ 



9A \ 



9 A 3 A 
18 1 9A 3A 



1 
18 



18 1 9A 
. A _|_ A / 

9A 18 T 9A / 



(69) 



Just as with graphs L 3 and K 3 , these two graphs have the 
same matrices because the Hamiltonians of the graphs 
KL 3 i and S4 commute; and, as we saw above, this pro- 
duces similar dynamics when we measure the walk in the 
corresponding final vertices v\ and v 2 . 

This is not the case, however, when the final vertex is 
v 3 or t>4 for these graphs. For S4, the *P and matri- 
ces for V3 and 1*4 can be obtained from those above by 
interchanging v\ with v 3 or W4. For KL31, however the 
matrices are 



<$ x (KL 3A ,v 3 ) 

/ 1 

10 

10 

\ 1 

Sj x (KL 3A ,v 3 ) 
( 



(70) 



V 





A J- A 

18 ' 9A 



X _ 1 
2 A 

A + 1 

"ilk 

A " r 2 





1 _ 3i 

A ± 2 

A 
j_ 

A 



A 



(71) 



The matrices ?$\{KL 3i \,Vi) and 9)\{KL 3 ^,V4) can be 
found by interchanging v 3 and U4 in the matrices above. 

The infinite hitting times for the L 3 and KL 3 .\ graphs 
can be understood to arise because the Hamiltonian of 
those graphs commutes with the Hamiltonian of a more 
symmetric graph. As we shall see in the next section, 
this fact can lead to infinite hitting times under certain 
circumstances. 



y x (KL 3 , 1 ,v 1 )=y x (S i ,v 1 ) 



(1 









1/ 



fi x {KL 3i x,vx) = fi\(Si,vi) 



( 



A 

_3_ _ 

4A 
J3_ _ 

4A 



-t + l ir 
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2 A 



4A 



A+8 



(66) 



A 4- i. 

4A ' 2 

A _ I + 

A + p 
A+8 



(67) 



V. INFINITE HITTING TIMES FOR GRAPHS 
WITH NON-CONNECTED COMPLEMENTARY 
GRAPH 

We will now look in a little more detail at infinite 
hitting times, which are one of the most surprising dif- 
ferences between classical random walks and quantum 
walks. In the classical case, if the graph is finite and 
connected, the probability to reach any vertex starting 
from any other is always 1. That is not the case for 
quantum walks. A sufficient condition for infinite hitting 
times was given in 37J: if the graph's Hamiltonian is 
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sufficiently degenerate, infinite hitting times will always 
exist. For continuous-time walks, any degeneracy at all 
is sufficient. (This is a sufficient but not a necessary 
condition because, as we've shown above, even graphs 
with nondegenerate Hamiltonians may have infinite hit- 
ting times.) We will now show that another sufficient 
condition for a continuous-time quantum walk to have 
an infinite hitting time is the non-connectedness of the 
complementary graph. Consider a graph T with n ver- 
tices and Hamiltonian Hr given by the usual expressions 
Eq. |T]) and Eq. (J5J). The complete graph K n with n ver- 
tices has the following Hamiltonian (in the basis spanned 
by the vertex states): 



1 



-1 



H K . = 



1 \ 

-1 



(72) 



\ -1 -1 ... 71-1/ 

This can be rewritten more succinctly as 

^„ = n(/-|Vo)(V'o|) = nPo ! 



(73) 



where |V>o) = ^ ELi \ k ) and p o = I - \ipo){i>o\- 

The complementary graph T c of a graph T is obtained 
by connecting vertices that are not connected in the orig- 
inal graph r, and removing the edges that are present in 
the original graph. Then it's easy to see that the Hamil- 
tonian of r c is 



H^c — H~k — H 



(74) 



Another observation is that the Hamiltonian of every 
graph r commutes with the Hamiltonian of the complete 
graph with the same number of vertices: 



[Hr, Hk„ 



0. 



(75) 



This follows from the observation that \ipo) is always an 
eigenvector of Hr with eigenvalue 0. As 



H r = P H r P , 



(76) 



([75]) is obvious. We can see from ((74|) that [Hr, H r -} = 0. 

Let us assume that the graph T is connected but the 
complementary graph r c is not, and consider the quan- 
tum walk on T c . Because T c is not connected, there are 
initial states that never reach a particular final vertex 
if the initial state includes only vertices which are not 
connected to the final vertex. Let's consider an initial 
state \ipi) that contains only vertex states that belong to 
one of the connected components of T c . Let us further 
assume that 



1 ™ 

(MM = -r y>l^> = o, 



which is always possible if the connected component has 
more than one vertex. Since is orthogonal to \iJjq) , it 



immediately follows that it is an eigenstate of Hk„ with 
eigenvalue n. 

If the final state \ipf) contains only vertices belonging 
to a different connected component of T c , the probability 
to ever reach the final state is 0: 



<V>/I' 



-itH T 



IVi) = ovt. 



(77) 



Putting all this together, we can see now that 

(l&/|e-* ffr M) = e -<*(^-*.) |^) 

= e- Un (^ f \e UH ^ \ifit) =0. (78) 

This proves that existence of infinite hitting time for the 
original graph T. 

We note that similar considerations may apply in some 
cases if the complete graph is replaced by a symmetric 
graph whose Hamiltonian commutes with with Hr- An 
example of this is the similarity of the dynamics of the 
KL31 and S4 graphs. 

Finally, we note that this sufficient condition for infi- 
nite hitting times is not particularly strong. For a graph 
with a large number of vertices, the complementary graph 
is almost always connected. 

VI. DISCUSSION 

We have examined continuous-time quantum walks 
and studied natural definitions for the hitting time. Af- 
ter considering different possibilities for introducing a 
measurement scheme, one of them emerges as a natu- 
ral one for the continuous case: measuring the presence 
or absence of the particle at the final vertex at Poisson- 
distributed random times, with an adjustable rate A. 
This is exactly equivalent to performing a particular type 
of weak measurement at frequent intervals, in the limit 
yielding continuous monitoring with time-resolution 1/A. 
Using this measurement scheme, we derived an analyti- 
cal formula for the hitting time which closely resembles 
the formula for the discrete-time case. 

This formula enables us to find a necessary and suffi- 
cient condition for the existence of quantum walks with 
infinite hitting times, namely that a certain superoper- 
ator pencil is not regular. In the case of finite hitting- 
times, the dependance of the hitting time on the rate of 
the measurement was studied, and the intuitive expec- 
tation for its behavior in the limits of weak and strong 
measurement rate was confirmed. In particular, as the 
measurement rate goes to infinity, the hitting time can 
diverge due to the Quantum Zcno effect. 

As in the discrete case, the symmetry of the graph 
plays a very strong role in the emergence of infinite hit- 
ting times. The graph symmetry group, if big enough, 
causes degeneracies in the eigenspectrum of the Hamil- 
tonian which in turn leads to the emergence of infinite 
hitting times for certain vertices. But this is not the 
only way in which symmetry can lead to infinite hitting 
times. Even when no degeneracy is present, symmetry 
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can cause some eigenvectors of the Hamiltonian to have 
zero overlap with some vertex states, as in the case of 
the L3 and KL3.1 graphs examined in section IV. This 
can be attributed to the fact that under the action of the 
group C2, the Hilbert space splits into symmetric and an- 
tisymmetric subspaces, and some eigenvectors from the 
antisymmetric subspace could have zero overlap with cer- 
tain vertex states. A further study exploring this idea is 
necessary to see if similar effects can occur for other sym- 
metry groups. 

Finally, in section V we show another condition for infi- 
nite hitting times. We have shown that the quantum walk 
on a connected graph can have infinite hitting times if the 
complementary graph is disconnected. This is in sharp 
contrast with the classical case, where every random walk 
on a connected graph will hit any vertex with probabil- 
ity 1 at long times. While this new condition is rather 
specific, it is possible that it can be generalized by re- 



placing the completely-connected graph with some other 
highly symmetric graph such that the Hamiltonian still 
commutes with the Hamiltonian of the original graph. It 
may be possible to explain any infinite hitting times on 
any graph in this way, giving a unifying view of the whole 
subject. It is clear that many questions remain, and that 
hitting times for continuous-time quantum walks are a 
very fruitful area of research. 
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